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ABSTRACT 


The equations of motion are developed and solved numerically 
for the trajectory of a spherical particle passing through a Prandil- 
tfeyer expansion fan. The effect of achangein y and n is shown, 
where y is the ratio of specific heats of the gas, and n/2 is the 


exponent in 
n/2 
tf 


2 (=) 
3 fe) 
the assumed viscosity temperature relationship. 

It is demonstrated that for particles the order of a micron in 
diameter, slip flow will exist, and a proposed correction to the drag 
to account for this discontinuous nature of the flow is investigated. 

The resulta are plotted showing particle trajectory profiles 


and the cormponents of the relative velocity or slippage velocity of 


the particle for flow deflection angles up to 70 degrees. 
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lL INTRODUCTIC?H 


Heterogeneous flow of a solid particle - gas rnixture has been 
a topic of increased interest in recent years. The addition of light 
metals to solid rocket fuels has focused attention on the non-equilibri- 
um nature of two-phase flows in expansion processes, and numerous 
papers have been written regarding this phenomenon. 

Early studies of two-phase flows were involved with investi- 
gations of water droplet trajectories around a variety of body shapes 
in connection with icing of aircraft surfaces. More recently, hcetero- 
geneous nozzle flow has been examined. Gilbert, Davis, and Altman 
(reference 1) related rocket thrust losses to solid particle sizes and 
solved linearized nozzle equations assuming Stokes drag law. MHliegel 
(reference 2) treated simultaneous velocity and thermal lag along with 
the effect of particle lag on gas properties in a one-dimensional noz- 
zie. Rannie (reference 3) showed that for most cases of interest, 
results of equal accuracy to numerical integration can be obtained 
through a perturbation analysis of temperature and velocity lags ina 
one-dimensional nozzle. Rannie also introduced slip flow corrections 
for drag and heat transfer. An excellent summary of recent advances 
in heterogeneous nozzle flow was presented at the American Rocket 
Society, Solid =ropellant Rocket Conference in January 1962, by R. 
i’, Hoslund (reference 4). 

Several problems relating to supersonic flow of two-phase 
systems have also been examined. Carrier (reference 5) investigated 
a solution to particle flow through normal shocks, and oblique shocks 


were treated by Morgenthaler (reference 6). 





The behavior of a single particle passing through a Prandtl- 


Meyer expansion furnishes an example of heterogeneous flow that can 
be analyzed in a straightforward manner, and serves to illustrate 
some of the problems and important parameters of two-phase flow in 
general. This paper presents the development of the equations of mo- 
tion of small spherical particles flowing with a gas that experiences a 
YFrandtl-Meyer expansion, and the numerical solution of these equa- 
tions. The effects of the temperature dependence of viscosity, changes 
in the ratio of ras specific heats, y, and slip flow correction to the 
drag coefficiont were investigated by varying these factors in succes- 
sive solutions. ‘fhe results of this investigation have direct bearing 


on the exhaust of an underexpanded rocket nozzle as the flow expands 





around the lip of the nozzle. 











Tl. DEVELOPMENT OF BASIC ECUATIONS 


An ideal two-dimensional Prandtl-Meyer expansion about a 
sharp corner in a semi-infinite channel is the assumed model for 
this investigation. A cylindrical coordinate system with origin at 
the vertex of the expansion fan is used to describe the position of a 
spherical particle passing through the fan. Oenoting the radial dis- 
tance from the origin to the center of the particle as , and the 


angular displacement of the particle from the initial characteristic 


* 





as (Fig. 1), the conservation of momentum equations can be 
written 
PP pa 2 coe o 
(mir) = mh od. Pr fF elou (us -r) (1) 
: Sony | Par 
ein = {Cv ew (ao . (2) 


It is assumed that the only force acting on the particle is a drag force 
which can be represented by Stokes drap law modified by the factor 

f . The radial and tangential components of the gas velocity are Ug 
and VU, respectively, and Be is the local value of the gas viscosity. 
The particle has mass m7 and radius oO. 

It is convenient to rewrite the equations with © the independ- 

ent variable. Choosing [= mr7?6 , the angular momentum, and 

Y- as the desired dependent variables, the time derivative can be 


written 


Equation 1 then becomes 
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Equation 2 is simply 
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Non-dimensional parameters are defined as 
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where = » & is the ratio of gas specific heats, and C 


is the theoretical maximum expansion velocity determined by the 


equation 
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The characteristic length, * » has the value 
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The subscript zero indicates initial conditions of the gas or particle 


as they enter the expansion fan. 


introducing equations 6 into equation 5, we obtain 


* 43 


s(e tz itty 6TT nk ‘ryt uy] = © 


and dividing by ae CITT M po 3 


Mage alt REL EYE + B]=0 


a“ Cc 





Equation 4 can be written in terms of the dimensionless parameters 


as 
d ‘ ‘4 OTe ol mer 
nC Vie 1 ORO SY Crees ') 
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The conservation equations of frictionless flow yield Ug Fy 
C 
ME} » and ae functions of ©. The equations may be write 
Cc ° 
ten in terms of our assumed coordinate system as 
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Equations 11 and 13 can be recombined to eliminate the term in the 
parentheses and show that 


Using this result in equation 14, V3 can be found as a function of 


(3 which can, in turn, be used in equation 12 to give 
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The value of the constant “& can be feund by evaluating equations 16 
at the initial conditions @=9% , Vy = GA. _, and using equation 


i4 to give 
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The quantity «can be found by applying equation 14 to show that 
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It is assumed that viscosity varies with temperature raised to a frace- 


tional power, hence 
ie 
Zz 


f- a ( 2 } (19) 
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or from equation 18 
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where |. will normally be between 1.0 and 2. 0. 


The equatiorsof particle motion (9 and 10) can now be written 





entirely in terms of the dimensionless variables R , “ , Oo | 





drag correction factor { , and the constants << and A. Letting 
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A value of [ for spheres in the subsonic slip flow regime has 
been proposed by \\. D. Rannie (ref. 3 ) which, for relative Mach num- 
bers less than unity and for Reynolds numbers less than about 0.5, 


takes the form 


i” ae 


= iste ae fs = i © (23) 
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where MM, and Ke are the relative Mach number and Reynolds 


number respectively. 
: Ms 
Ke 


path to the sphere diameter,can be evaluated as follows: 


» which is related to the ratio of molecular mean free 
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From equations 15 and 16 
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Isentropic relationship between density and termperature along with 


equation 19 gives 
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Using equations 20, 25, and 26 in equation 24, we get 
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and equation 23 becomes 
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The drag correction factor cannot be written in terms of the dimen- 
sionless similarity parameters Ro and A , but is dependent on the 
particle size and initial conditions. 

It is convenient to relate the trajectory cf the particle to that 
of a particle having zero relative velocity, i.e., moving along a gas 


streamline. Desionatines this oath as a Se0lr eC we can write 
o “> p 3 3 ) 


el ay. meee 5 
dt 1“ de dt 


. > - 

















a8 
¥ 
=, 
- Je 





AY fm 


a 





and since Go - 2 we have 
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Integrating this cquation we obtain 
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From the definitions of AR and A it follows that 


Va Fs A, 
nue Fe, No 





which, when combined with equations 16 and 30, gives 


{ 

4 -) 
Ag _ [cos & S 
Ao  \COS(AOE+ &) 


Two new dimensionless parameters are defined 
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from which it follows that 
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| Also from equations 6 and 16 
Roe cone (38) 
Writing equations 21 and 22 in terms of X , Y parameters 


using equations 33 through 355, we obtain the result 
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It is assumed that the particle is moving with no relative velocity prio 


to reaching the first characteristic line of the expansion. The initia! 


conditions are then 
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Equations 28, 39, 40, and 4] constitute the complete statement of 


the problem. 
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Ill. AUXILIARY ECUATIONS 


In order to specify the trajectory, the velocity as well as the 





position of the particle should be known. The important parameters 


in terms of X and Y are then 
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which, by 33, 34, and 16, is 
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The relative velocity or slip velocity of the particle relative 


to the absolute gas velocity is simply 
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The relative Mach number is then 
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IV. SCLUTION OF ECUATIONS 


Since equations 39 and 40 are non-linear, a numerical solu- 
tion is indicated for arbitrary values of RR, . Foreither FR, 77 !.0 
or R,<< 1.0. , the possibility of a small perturbation solution 
exists; however, this possibility was not investigated. 

Equations 21 and 22 were first solved utilizing a desk calcu- 


lator and a simple iteration procedure where 


d ( di (47) 
Gal= (e+ (Rl ee 
Ae 2[ “ay ee, eaal®? - 


For these calculations it was assumed that {=)].0 (Stokes drag law 
applies) w-!.4, mM=2.0 , and Mo-— /2’ . Various values 
of RK, between 2.1 and 100 were investigated; however, it was 
found that the iteration interval, A 6, hadtobe decreasedas AY 
increased to keep the process convergent. Tor AY, very small, the 
iteration interval hac to be made very small again to get meaningful 


result (proper divergence frorn the zero deflection line). 
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To eliminate much of this difficulty, equations 39 and 40 were 
developed. A similar iteration procedure to that of equations 47 was 
initially used to investigate the new equations, and it was found that 
an iteration interval of 0.901 radians gave satisfactory results for FR, 
between 0. 01 and 190. 

Equations 39 and 40 appeared to be the more satisfactory of 
the two sets of equations for programming on a disital coraputer, since 
the starting solutions of numerical methods are closely related to the 
Sirnple iteration procedure. 

A Fortran statement for the solution of equations 39 and 40 
was written utilizing a subroutine prepared by Mr. K. H. Redcner of 
Computer Sciences Corporation. The subroutine applies the method 
of Runcge- Xutta-Gills as a starting solution, and the Adams-Moultoa 
predictor - corrector formulas as a continuing solution. The facilities 
of the Computer Center at the California Institute of Technology were 
used in writing the program, and the IRM 7090 computer at the Jet 
Fropulsion Laboratory was used for the actual solution. A Fortran 
statement is included as Appendix A. 

To investigate the general effect of changes in ¥Y and hn , 
Stokes drag law was assumed to apply, and solutions were obtained 
for six values of ARQ, between 9.01 and 50 for three sets of assumed 


constants. 
1) M > = /2 ; wie 1.4 5 ce O 
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The exhaust conditions of a typical large rocket were assun 


as initial conditions for the investigation of the effect of allowing 







to vary according to equation 25. The initial conditions chosen were 





We (Ona 

Ho =~ 9.55 X10" eae 
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me) = le oe 

A ee he. 


Solutions were obtained for A, of 0.1, 1.0, and 50 for two par- 


ticle sizes, © = 4 micron and ¢o = 2.0 microns. 
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V. RESULTS AND DISCUSSION 


The results of this investigation are plotted on Figs. 3 through 
10 for 0° to 70° of flow deflection ancle. Flow deflection angle, a, 
was chosen as a more meaningful parameter than ¢ for plotting and 
analyzing the results since it is determined solely by the physical 
boundaries of the flow. The fan angle and deflection angle can be re- 
lated by the usual Prandtl-Meyer functions. 

Figure 2 is a graphical representation of equation 6 giving Fy, 
as afunction of Yo and C for particles having a density of 
210 lb/ft” (as for ALC, ). Initial gas characteristics plotted are for 
air at 100°F and 1000°F and for the assumed rocket exhaust conditions 
given in the previous section. AX, is not greatly affected by a change 
in initial gas conditions, but is primarily a function of fm and CT . 
For particles between 1.0 u and 4.0 41 in diameter and initially dis- 
placed by several cm. from the flow boundary (as in the case of a 
rocket exhaust), A, will normally be greater than 10. 

The data plotted in Figs. 3, 4, and 5is based on the assump- 
tions that Me=/2. » ozl.t , m=z 2.0, and Stokes drag law ap- 
plies with the factor {= /.O . For these conditions, particles with 
R,.?50O- will have little deviation from the streamline, and slip 
velocities will be low. Sor AR. 4 0.0] there will be large tan- 
gential slip and very little deviation frorn its original flow direction 
through the first 25° of gas deflection. Radial slip is small compared 
to tangential slip, being a maximum of .164 for AR, =0.0/ at about 


25° flow deflection. Near this point the maximurn relative Mach 
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number is .993 and maximum Reynolds number for a 4. 0 y particle 
is 858. For R,=10, the maximum radial slip is less than one per 
cent while the maximum Mach number is only 0.151] and maximum 
Reynolds number for 4.9 y particle is 2.U2. ‘fhe Mach numbers and 
Reynolds number both drop sharply on either side of the maximum 
values. 

The assumed viscosity - temperature function, equation 19 

n 

ae ( = ) 2. (19) 

fo To 
is a well-established approximation that applies over a limited tem- 
perature range when / is assigned the constant value appropriate to 
this range. By allowing to take its expected limiting values, its 
effect on the solution can be observed. Figures 6 and 7 show curves 
for which Stokes drag, Y=/.4¢ , and M, = V2 where assumed for 
N=l.0 and h=2.0o . A small change in trajectory path and tan- 
gential slip velocity and a negligible change in radial slip velocity 
were observed. Thus, little error is introduced by the selection of 
an average value of /) over an extended ternperature range. 

Achange in Y from 1.4 to 1. 2 does have a fairly large ef- 
fect on the trajectory path which is due principally to the change in 
the gas streamline. The tangential slip velocity is increased also, 
which tends to increase the path deviation from the streamline, but 
this is a small effect. Radial slip velocity is not appreciably affected. 

The drag correction factor, | , which allows a departure from 
Stokes drag law to compensate for slip flow or free molecular flow 


effects, was proposed by Rannie (ref. 3) as an analogue to a similar 





Me 
development for heat transfer coefficients. Although it has not been 
verified experimentally, it does give proper results for the extreme 
values of Ms yl O and “%s ¢¢ ;.0 » and should give a good 
Fe e R = M 
approximation for intermediate values. The factor —° canbe 


Re 
shown to be (ref. 7 ) 


where | is the mean free molecular path of the pas. Thus, if 


Ms > > (fee , free molecular flow exists, andif M/S cc to 
Re Re 


continuum flow exists. The area between these limits is the slip 
flow regime. 


From equation 27 we have, for S=l-2 and n= /7em 


pa. K ( ea ) Ag 
Re Cos(pO+aX)/ | (oa) 


For air or rocket exhaust gases at moderate temperatures, and for 
particles the order of a micron, kK will lic between 0.03 and 0.3. 
Thus, a particle of this general 3ize, passing through a Frandtl- 
Meyer expansion, will start in the slip Now regime and pass into free 
molecular flow at large expansion angles as the cnclosed cosine func- 
tion in ecauation <8 becomes large compared to unity. 

Figures 5, 9, and 10 show the effect of f{ on the trajectory 
paths and slin velocities for 1.2 micron and 4.9 micron diameter 
particles. The deviation from the uncorvected Stokes drag law trae 


jectories is large for particles as small as 4.9 micron diameter, 

















-20- 
and increases markedly as the particle size is reduced to 1. 0 micron. 
In the case of a 1.6 micron diameter particle, an R,.=50 does no 
closely follow the streamline as it did for Stokes drag, and Fr = © 
now shows little deflection from its oripinal flow direction. For 
R.,=0O.|1 a maximum relative Mach number of 0.924 and maximum 
Reynolds nurnver of 6.73 is reached for ea 4.0 micron particle. This 
compares to maximums of 0.853 and 6.51 respectively for Stokes 
drag. 

The maximum values of Reynolds numbers obtained are greate 


than the limits set on equation 23 for the determination of £ . A 


Reynolds number of & 0 would change the equation to 


M*s 
— = ©. 6 Ou mawo > VY — 3 
f i 


(ref. 3}, which would make a 23 per cent errorin { for a Maci num- 
ber of 0.95. While this is an appreciable crror, itis effective over 
a small portion of the total expansion and diminishes rapidly on both 
sides of maximum A, , so it should not greatly affect the overall 


results. 
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VI. CONCLUDING REMARKS 


In the flow repime where Stokes drag law applies, the tra- 
jectories of spherical particles passing through a P'randtl- Meyer ex- 
pansion fan can be specified by a single dimensionless parameter, Re ’ 
which is a function of initial gas conditions, initial radius, and particle 
diameter. If the Reynolds number is greater than about 0.5, the tra- 
jectories are also a function of Re. When the mean free path of the | 
gas molecules is of the same order or large compared to the particle 
diameter, Stokes drag law must be modified to correct for the discon- 
tinuous nature of the flow, and the trajectories then become functions 
of three parameters, Ro , Réweand Ma me. 

It was found that the temperature dependence of viscosity need 
not be precisely known to ensure satisfactory results, but changes in 
the specific heat ratio, y , do appreciably :ffect the trajectory paths. 

The radial velocity lags are in general very small compared 
with the tangentiai velocity lags, and can, in many cases, be neplected 


completely. 
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APPENDIX A 
FORTRAN STATEMENT 


CIMENSICN Y(1C), YOCT(3), Cl(1C), C10) 


CIMENSILON Y€10), YCCT(3), C(lC), C10) 
5CG FCRPAT (7F10.4) 
51C FCRMAT (1CK1 GAMMA = El2.6, 13 MZERC = E136 


19F N = F522, 13H RZERC = E13.6, 13H DELTA = E13.6//) 


1C READ INPUT TAPE 5,5CO, CGAPMFMA, CMZERO, CN, RZERO,DELT, CMRE, 
WRITE CUTPUT TAPE 6, 510+ CGAPMMA, CMZEROs CNye RZEROe DELT 


Y(l) = 1.€ 
Y(2) = 1.0 
Y(3) = 0.€ 
T = 0.0 
KAY = 1 


CLAMCA = SCRTFI(CGAPMMA-1.0)/(CCAPMMA4+1.0)) 

B = (1.0/(CLAMCA##2) )-CNH1.C 

ARG = SOQRTF (1.-0+CLAMDA#CLAMDA#(CMZERC#CMZERC—-1.0)) 
ALPHA = ACCSF(1.0/ARG) 


Ctl) = (1-C/CLAMECA)-CLAPDA 

Ct2) = C(1l)/CLAVCA 

C(2) = CCSF (ALPFA) 

C(4) = CLAVMCA#(C(3)##(B-1.0))*#RZERO 
C(5) = 2e-U/CLAPDA 

Clteye- CRRE / CC 


hRITE CUTPUT TAPE 6,5CO,CLAMCA,B,ARG,ALPHA 

WRITE CUTPLT TAPE 655CCoCE1) C12) 9013) C14) C15)5 Cl6) 
CALL CEC(K,3e7TsYeYDGT,.C1,5.CE-7) 

GC TC (100,20C,2C0,30C),K 


1CC Cll) = CLAMCA#T+AL PRA 
Ci2) = TANFI(C(1)) 
EiSo7e= CESFIC(1)) 
C(4) = (1.C/(0(3)##B)) 
CIS) = 1.C/(1.0+C(6) * ((C(3)/0(3)) ##(8-1.0))) 


VEETUL = vyi3) 

VECT C29 5= (C(5)) #v¥(3) sD(2) — Cl2)4Y( 1)*Di(2)*D{(2Z) 

I-(CU4) #V(3)/00VOL) ##3) #¥(2)8#4¥(2)))*#D(4) * DIS) 

YDCT(2) = -CLlL)#V(2)#C(2)4C(4)#D(4) (10/00 V01L) ##3) #¥(2)) 

P=1 607 (¥il)*Y{(1))) #C(5) 

IF (KAY-2) 150, 160, 16C 
160 CALL CEQ2 
15C wRITE CUTPUT TAPE 6, 500, Cl), Dt(2), C(3), D(4),D(5) 

WRITE CUIPLT TAPE 6, 60C 


6CC FCRMAT (//6F THETA,9H X~, L4H Yo 
114k Z, 16H V/VS,14H U/US, 
214+ R/RC, 13h MS, 13H RE//) 
KAY = 2 
GG TC 160 

2cc ¥(4) = Y(1)#Y(2) 
Y(5) = Y(4)#(1-0-CLAMDA#YDCT(1L)/(Y¥Y(1)#D(2))) 
Y(@é) = ((C(3)/003))##(C(2)4+1.0))/Y (1) 
Y¥(S) = ((1-C-Y¥(4)) ##2)4(°0(01-C-¥(5))#D(2)/CLAMDA) # #2) 
Y(7) = SCRTF(Y(9)) 
Y¥(8) = Y(7)#((C0(3)/C(3))*##(B-1.0)) 


wRITE CUTPUT TAPE 6¢ 610, Te (Y(JACK),JACK = 1,8) 
610 FCRMAT(/F7.3,8E14.6) 
205 ITFLT—3.0)2160,210,220 
21C CALL CEQL 
GEST 205 
220 GC TC 10 


CIMENSICN YO(1C), YDCT(3), C10), C10) 


3CC WRITE CUTPLUT TAPE 6% 62C 

620 FCRMAI (22K ERRCR RETURA FRCP CEQ) 
GC TC 10 
ENC (1 s0e0e0e0r0eCe0e0elsOeCe0eCe0) 























